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CALCULATION  OF  PENETRATING  RADIATION  OF  THE  AIR  DEFENSE  STRUCTURE 
Kang  Ning 

The  fundamental  factors  that  kill  and  injure  people  in  a  nuclear  detonation 
are  shock  wave,  photo-radiation  and  nuclear  radiation.  During  a  detonation  of 
a  nuclear  device  at  a  height  lower  than  30  kilometers  in  the  atmosphere,  the 
energy  of  photo-radiation  occupies  35  percent  of  the  detonation;  however,  50 
percent  of  the  nuclear  detonation  is  expended  in  forming  a  shock  wave.  The 
remaining  15  percent  of  the  energy  is  released  in  the  form  of  nuclear  radiation. 
Although  the  energy  of  nuclear  radiation  occupies  a  small  portion  of  the  total 
energy  of  detonation,  the  radiation  has  very  high  penetrating  power.  Therefore, 
at  a  relatively  great  distance  from  the  detonation  center,  even  the  killing  and/or 
injuring  function  may  exist  inside  a  shelter.  Therefore,  the  study  of  protection 
from  nuclear  detonation  is  a  difficult  task  for  technical  personnel  in  underground 
construction. 

The  calculation  of  underground  construction  providing  protection  from  nuclear 
radiation  includes  three  parts;  the  top  cover,  wall  and  hole  opening.  It  is 
required  that  the  incident  dosage  does  not  exceed  the  allowable  dosage  of  50 
roentgens.  However,  in  the  present  papers,  there  is  a  lack  of  clarified 
description  on  this  problem;  in  addition,  points  in  the  calculations  still  have 
room  for  further  discussion.  Mostly,  only  incidence  of  nuclear  radiation  along 
the  normal  direction  is  considered,  without  considering  oblique-direction  nuclear 
radiation.  Although  in  some  cases  oblique  incidence  is  considered,  oblique 
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incidence  after  the  radioactive  cloud  rises  is  not  considered.  For  the  afore¬ 
mentioned  situations,  this  paper  presents  the  problem  of  protection  from 
nuclear  radiation  for  underground  construction.  Readers*  corrections  are 
welcomed. 

I.  Characteristics  of  Nuclear  Radiation  (Penetrating  Radiation)  and  Decay  of 
Radiation 

1.  General  characteristics  of  penetrating  radiation 

Nuclear  radiation  is  composed  of  alpha  particles,  beta  particles,  gamma  rays 
and  three  types  of  neutron  streams.  The  effective  range  of  alpha  particles  in 
the  atmosphere  is  only  several  centimeters;  the  effective  range  of  beta  particles 
(high-speed  electrons)  is  still  only  several  meters  in  the  atmosphere.  The  range 
of  gamma  rays  (high-energy  electromagnetic  radiation,  similar  to  X-rays)  is  as 
far  as  thousands  of  meters  in  the  atmosphere;  the  range  of  neutrons  (particles 
released  after  nuclear  fission  and/or  fusion)  is  hundreds  of  meters.  It  is 
apparent  that  the  calculation  of  penetrating  radiation  is  mainly  the  study  of 
ganana  rays  and  calculation  of  the  neutron  stream. 

However,  in  actual  structural  calculations,  we  do  not  consider  all  gamma  rays 
and  neutron  streams  released  following  the  detonation  of  a  nuclear  weapon.  For 
detonation  of  fission  bombs,  according  to  the  emission  time  the  gamma  radia¬ 
tion  is  divided  into  instantaneous  gamma  radiation,  short  lived  gamma  radiation, 
and  gamma  radiation  lasting  a  number  of  seconds.  As  revealed  in  studies,  the 
major  radiation  dosage  in  the  detonation  area  is  radiation  lasting  a  number  of 
seconds,  and  next  important  is  the  short  lived  radiation.  The  longest  acting 
time  of  these  radiations  on  the  ground  target  is  15  to  20  seconds.  Neutrons 
from  nuclear  fission  can  be  divided  into  fast  neutrons,  medium-speed  neutrons 
and  slow  neutrons.  Since  the  dosage  of  slow  neutrons  is  smaller  than  that  of 
fast  neutrons,  and  at  the  same  time  neutrons  are  released  in  a  relatively  short 
period  of  time,  the  casualty  function  of  neutrons  is  generally  caused  by  fast 
neutrons.  It  is  apparent  that  the  calculation  for  a  protective  structure  only 
considers  gamma  radiation  lasting  a  number  of  seconds,  short  lived  gamma  radia¬ 
tion  and  fast  neutrons.  However,  the  penetrating  radiation  dosage  released  by 


a  fusion  bomb  (thermonuclear  bomb)  includes  some  ultra  high-speed  neutrons 
beside  gamma  rays  and  high-speed  neutrons. 


2.  Decay  process  of  penetrating  radiation  (refer  to  two  flow  charts  in  the 
following) 
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Key:  (1)  Nuclear  fission;  (2)  Fission  fragments;  (3)  Fast 
neutrons;  (4)  Gamma  rays;  (5)  Slow  neutrons;  (6)  Capture; 

(7)  Excited  nucleus;  (8)  Radioactive  nucleus;  (9)  Fission 
gamma  rays;  (10)  Captured  gamma  rays;  (11)  Gamma  rays  with 
inelastic  scattering;  (12)  Instantaneous  fission  gamma  rays; 
(13)  Residual  radiation;  (14)  Absorption;  (15)  Photo-excited 
neutrons;  (16)  Electron-pair  effect;  (17)  Positrons  and 
[Key  continued  on  following  page]  , 


[Continuation  of  Key  from  preceding  page] 

electrons;  (18)  Cascade  radiation;  (19)  Photoelectric  effect; 
(20)  Photoelectrons;  (21)  Compton-Wu  (Youxun)  effect;  (22) 
Nuclear  scattering  gamma  radiation;  (23)  Nuclear  released 
electrons;  (24)  Penetrating  radiation  entering  underground 
works  (air-defense  structure). 


III 


1 

(3) 

*  *  * 

(4) 

He*Qa|&?  | 

{* 

I* 

m. 

;£ 

(8)  * 

it* 


(5) 


*  *  * 


R  (6 


(7)  |  ***** 


|C11) 


(9)| 

T 

ft 

t* 


Kim* 


(1C) 


I  & 


It 

tft 


8\JfeTi£  ISHMIt  (22) 


Key:  (1)  Nuclear  fusion;  (2)  Ultra  high-speed  neutrons; 

(3)  Fast  neutrons;  (4)  He4  (alpha  particles);  (5)  Slow 
neutrons;  (6)  Capture;  (7)  Excited  nucleus;  (8)  Gamma  rays 
with  inelastic  scattering;  (9)  Captured  gamma  rays;  (10) 
Decayed  gamma  rays;  (11)  Radioactive  nucleus;  (12)  Absorption 
(13)  Electron-pair  effect;  (14)  Positrons  and  electrons; 

(15)  Cascade  radiation;  (16)  Photoelectric  effect;  (17)  Photo 
electrons;  (18)  Compton-Wu  (Youxun)  effect;  (19)  Nuclear 
[Key  continued  on  following  page] 
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[Continuation  of  Key  from  preceding  page] 

scattering  gamma  rays;  (20)  Nuclear  released  electrons;  (21) 

Residual  radiation;  (22)  Penetrating  radiation  entering 
underground  works  (air-defense  structure). 

II.  Calculation  of  Protection  by  Underground  Works  From  Gamma  Rays 

1.  Calculation  of  protective  layer  (including  cover  and  wall) 

When  the  gamma  rays  transmitted  from  the  detonation  center  to  underground 
works  pass  through  the  protective  layer  to  enter  the  underground  works,  the 
intensity  of  the  gamma  rays  decays  continuously  because  of  mutual  actions  among 
atoms  (in  the  forms  of  photoelectric  absorption,  Compton-Wu  (Yunxun)  scattering 
and  electron-positron  pair) .  Refer  to  YUANZI  BAOZHA  DI  GAMMA  FUSHE  [GAMMA 
RADIATION  FROM  ATOMIC  DETONATION],  page  33.  When  gamma  rays  have  a  vertical 
incidence,  the  calculation  formula  is 

Dr-DaKua rDne-‘  (j) 

In  the  equation,  D^  is  the  dosage  (in  roentgens)  of  gamma  rays  after  passing 
through  the  protective  layer;  D  is  the  dosage  (in  roentgens)  of  incident  gamma 
rays;  e  is  the  base  of  the  natural  logarithm;  y  is  the  total  decay  coefficient 
(obtained  from  Table  1)  of  gamma  radiation  exerted  by  the  protective  layer;  ay 
is  the  decay  coefficient  of  the  softest  scattering  of  the  protective  layer 
(the  value  of  the  decay  coefficient  can  be  obtained  from  Fig.  1);  K  is  the  decay 
coefficient  of  the  dosage  of  gamma  radiation  exerted  by  the  protective  layer;  and 
x  is  the  thickness  (in  cm)  of  the  protective  layer. 

In  the  situation  of  oblique  incidence,  decay  of  gamma  ray  intensity  is 
faster  than  in  the  normal  direction  (vertical)  incidence  through  the  protective 
layer,  because  the  passage  of  the  gamma  rays  is  longer;  in  other  words,  the 
thickness  of  the  protective  layer  along  the  direction  of  the  ray  beam  is 
relatively  increased.  From  0.  I.  Lieyisiji  [transliteration  of  Russian  name], 
at  that  time  the  dosage  decay  coefficient  KY  can  be  expressed  by  the  following 
equation: 
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Table  1. 


[a)_  ,  (b)  *«****MM*-'> 


« 

w*m  : 

(d) 

*** 

0.033^  1 

0.027  (-Cj 

(e) 

±  - 

0.08  | 

0.05 

(f)  mm± 

|  0.11 

0.07 

Cg) 

* 

0.33 

0.25 

(h) 

« 

0.80 

0.39 

m  m(n  «,-*?»!  Cj) 


Key:  fa)  Material;  (b)  Effective  decay  coefficient 
y  (cm-*);  (c)  Electron  megavolts;  (d)  Wood;  (e)  Earth; 

(f)  Concrete;  (g)  Iron;  (h)  Lead;  (i)  Remarks;  (j) 

Photon  energy. 

*»<*.*.«•>»  a<f.e#,Jt)Kr<lfe,)  (*> 

In  the  equation,  K^(2,  eQ)  is  the  dosage  decay  coefficient  of  gamma  rays  (with 

e  as  the  photon  energy)  while  passing  through  a  protective  layer  (with  thickness 

®  -Ux  -x/d 

l)  along  the  normal  direction;  K  is  the  calculated  as  K=e’  =2  '  ;  and  a(<l>,eo,K) 

is  the  coefficient  of  oblique  incidence  (refer  to  Table  2)  when  the  incident 

angle  is  between  0°  and  30°.  a($,  e  )*1;  i.e.,  the  oblique  thickness  of  the 

protective  layer  can  be  considered  the  same  as  the  vertical  thickness. 

As  is  well  known,  the  radioactive  cloud  (major  radiation  source  of  gamma 
radiation)  formed  after  a  nuclear  detonation  rises  with  a  definite  speed  (vc*72 
to  100  m/sec).  The  radiation  dosage  of  gamma  rays  acting  on  a  protective  layer 
of  fortification  and  the  minimum  thickness  of  (radiation  passing  through)  the 
protective  layer  vary  continuously.  For  different  oblique  angles,  there  is  a 
different  function  of  the  protective  layer.  The  smaller  the  oblique  angle,  the 
greater  is  the  relative  thickness  of  the  covering,  and  the  thinner  the  relative 
thickness  of  the  wall.  Conversely,  the  greater  the  oblique  angle,  the  thinner 
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is  the  relative  thickness  of  the  covering  and  the  greater  the  relative  thickness 
of  the  wall  (Fig.  2). 
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Table  2.  Coefficient  aj($,  eQ,  K)  of  oblique  incidence. 
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1.41 
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2.25 
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1.00 

1.00 

1.00 

1.00 

1.00 

1.05 

1.20 
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1.01 

1.23 

1.40 
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2.05 

2.50 

1  1.03 
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2.00 

2.50 

3.50 

4.10 

5.60 

1  1  1.46 

2.38 

2.00 

4.80 

6.30 

8.00 
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1  !  1.30 

3.40 

4.30 

8.30 
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2.9o 
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1  9.20 
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132 

•  1  !  4.60 
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43.8 
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132 

sir 

1  !  3.30 

;  1 

17.4 

24,0 

78.0 

211 

770 

!  1  {  10.3 

28.7 
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132 

205 
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1880 

j  1  i  13.8 

41.0 

57.0 

23.0 

382 

1500 

4800 

!  1  ;  24.4 

75.0 

88 

402 

342 

3360 

11200 

- - 

•  1 

38.7 

107 

136 

704 

1140 

8050 

27200 

I  1  35.7 

150 

210 

1240 

20400 

23000 

65800 

>.  Nim~co5T 

(b)  -*nnn.xm*m'*wf&rn 


Key:  (a)  Remarks;  (b)  This  table  is  based  on  Table  21  of  YUANZI  BAOZHA 
GAMMA  FUSHE  [GAJWA  RADIATION  OF  ATOMIC  DETONATION],  and  "uses"  a  hori¬ 
zontal  protective  layer;  (c)  In  the  calculation  of  dosage  decay  of 
penetrating  radiation,  $  is  calculated  by  extrapolation. 

At  that  time,  the  dosage  of  gamma  radiation  decayed  after  passing  through 
the  protecti'  j  layer  can  be  calculated  according  to  the  following  equation: 
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Fig.  1.  Variation  of  a  of  earth  Fig.  2.  ai  (oblique  angle);  $  (incident 

and  concrete  with  R/H.  angle). 

Key:  (1)  Rising  direction  of  radioactive 
cloud;  (2)  m/sec. 


Df*  £  <  »  > 


Expressed  with  half-value  thickness: 

d,- 


(  4  > 


If  we  simplify  Eq.  (4)  into  a  form  similar  to  incidence  along  the  normal  direction: 

Dr  =  afD„ a»'<  (4') 


In  the  equation,  d  is  the  half-value  thickness  of  the  material  (Table  3);  and  b 
is  the  decay  coefficient  of  oblique  incidence.  .  «  • 


Tables  4  and  5  show  the  oblique-direction  decay  coefficient  of  concrete  (reinforced 
concrete)  to  an  atomic  bomb;  x^  is  the  oblique  thickness  of  the  protective  layer 


in  different  time  intervals. 
For  cover: 


For  wall : 


,  .  H. 

-  (5> 

’  co^, 

»• 

K»c  rv'(H-»t,ll)*  +  R*  . 

'  rtn*. 

~  R' 
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In  the  equations,  Hn  is  the  thickness  (in  meters)  of  the  cover  protective  layer; 

is  the  thickness  (in  meters)  of  the  wall  protective  layer;  H  is  the  altitude 
(in  meters)  of  the  atomic  detonation;  R'  is  the  distance  (in  meters)  from  the 
underground  works  to  the  projecting  point  of  the  detonation  center;  j--t^, 

t^  i  is  the  intersection  angle  between  the  instantaneous  incident  gamma  rays  and 
the  normal  line  of  the  protective  layer;  and  t^,t^  ^  indicate  time.  The  action 
time  on  the  underground  works  by  gamma  rays  is  divided  into  n  sectors.  t^(t^  j) 
is  the  time  duration  between  the  instant  at  the  end  of  the  i(i-l)-th  sector  and 
the  instant  that  the  atomic  bomb  begins  to  detonate.  This  time  duration  is 
approximately  15  seconds.  Refer  to  Table  6  for  dose  percentages  released  at 
various  instants. 


Table  3.  Half-value  thickness  of 
materials. 


Table  4.  Oblique-direction  decay 
coefficient  b  with  cover. 


(a)  tmStt 

(d)  * 

(e) 

(f)  * 

(h)Rfl:?P5,  fcsm 

(i)»jmg, 

(jj  tews,  k 

(1)««± 

(m) 


<£'!§*)  i  d. 

>)/  ~tcr50 


0.5  1.0  1.5  2.0 

0,710  ^  0.390 | 0.220  0.14 

0.500  !  0.250 ’ 0.090  0.065 

j - j - ; - 

1.460, 0.200  I  0.061  0.052 


Key:  (a)  Material;  (b)  Weight  of  unit  f 

volume  (grams/cubic  centimeter);  (c) 

Half-value  weakened  thickness;  (d)  Snow; 

(e)  Wood;  (f)  Ice;  (g)  Soil,  sand;  (h) 
Weathered  gravel,  soft  shale;  (i)  Sandy 
shale,  schistose  sandstone;  (j)  Granite, 
gneiss;  (k)  Dressed  masonry;  (1)  Con¬ 
crete;  (m)  Reinforced  concrete;  (n) 

Iron. 


Note:  For  Tables  3  and  4,  data  are 
calculated  for  q  equal  to  8000  tons, 
and  H  equal  to  400  meters  altitude 
with  atmospheric  detonation  of  the 
atomic  bomb.  If  q  is  greater  than 
[sic]  8000  tons,  application  of  data 
in  both  tables  provides  higher  factor 


Table  5.  Oblique-direction  decay  coefficient 
b  with  wall. 


i  0.5 

!  o.tr  1 

1  1.0 

2.0 

HCr/do 
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0.48  i 

!  0.39  { 

0.17 

0.030 

10 

0.40 

|  0.30  1 

0.13 

0.01S 

12 

0.33 

0.25 

0.093 

0.0067 

IS 

0.30 

0.19  ' 

| 

0.61 

i 

0.0036 

18.5 

0.19 

:  o.i6 

! 

0.038 

0.00071 

20 

!  0.18 

;  0.12 

0.033 

'  0.00056 

Table  6.  Radiation  dosages  released  at  various  instants  with  atmospheric 
detonation  of  an  atomic  bomb. 
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20 

i  30 

150 
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5000 
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i 

|  30 

i 

:  iso 

i 

5000 

0—20 

1  o  ! 

i 

0 

1 

0 

i  • 

0.4 

1  0.03 

0.04 

j  1.35  ' 

2.2 

20—40 

;  0.025  j 

0.03 

j  0.40 

0.35 

2.2 

0.13 

0.16 

0.28 

|  1.2  : 

3.8 

40—60 

'  0.13  i 

0.16 

|l.28 

1.2  j 

3.8 

0.55  j 

|  0.7  1 

1.0 

!  30  1 

50 

60—80 

1  0.55  ; 

0.7 

1.0 

3.0 

5.0 

3.0 

...i 

3.0 

I 

I  5.0 

7.5 

80—100 

!  2.0  | 

2.3 

» 

5 

7.5 

IS 

15  ' 

15 

i  is 

20 

Key:  (a)  Types  of  protective  layer;  (b)  For  wall;  (c)  For  cover;  (d)  Mean 
calculated  time  (in  seconds)  of  the  assumed  percentage  of  radiation;  (e)  TNT 
equivalent  (in  kilotons);  (f)  Radiation  percentage. 

III.  Calculation  Examples  of  Gamma  Ray  Defense  Provided  by  Air  Defense  Structure 

There  is  an  annex  type  air  defense  structure  (six  stories  with  3.2  meters 
as  the  height  of  a  story);  the  TNT  equivalent  of  an  atomic  bomb  is  8000  tons 
with  a  detonation  altitude  of  400  meters.  The  detonation  occurs  at  the  lateral 
rear  side  of  the  structure.  The  distance  between  the  structure  and  the  projectng 
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point  R*  of  the  detonation  center  is  400  meters.  Calculation  is  required  to 
check  the  structure's  capability  of  defense  against  gamma  radiation. 


1.  Consideration  is  given  to  the  incident  dosage  of  gamma  radiation  on  the 
air  detonation  effect: 

2 

The  calculation  formula:  0  *AK/R  (7) 

yo 

In  the  equation,  K  is  the  relationship  coefficient  of  TNT  equivalent  q;  and  A  is 
the  influence  coefficient  of  photon  energy  of  gamma  rays,  considering  the  air 
detonation  density  p. 

From  the  dosage  distribution  (Fig.  3)  of  gamma  radiation  during  atomic  detona 
tion  of  a  Hiroshima  bomb,  we  can  obtain  the  following  simplified  formula: 

When  R<R. ,  A*l;  when  R.^R^R_;  A=e~R^25®;  as  recommended  in  some  manuscripts, 

...  .  -R/300 

then  A=e 


From  the  physical  concept  of  studying  the  function  of  the  shock  wave,  0.  I. 
Lenpunskiy  derived  the  following  equation: 

For  detonation  in  the  atmosphere:  K  •  t.**!*^**2*^®  •••***» *'«>  (8) 

Simplifying  Eq.  (8),  we  obtain  ..<•*) 

Then  the  value  of  K  can  be  obtained  from  Fig.  4. 

For  ground  detonation:  K=2.8xl0^q(l+0.3q®*65)  (9) 

When  R>R2,  a* 

K-S.SxlOVrCl  *0.0$<aq)t ■*]  UO) 

Or,  we  simplify  it  into  K=5,5xl0^q  (10*) 

For  ground  detonation,  a=2;  for  air  detonation,  a=l.  R  is  the  distance  from  the 
air  defense  structure  to  the  detonation  center. 

r=/h2+r’2  =  / 4002+4002* 565  meters. 

This  topic  is  ltl<R<IL  A  ««-*'*• 


From  Fig.  4,  we  obtain  K=1.8xl0  .  The  incident  dosage 
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Dy0* [ (1 . 8x10  ) /565  ] xO .208=10,700  roentgens . 
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Fig.  3.  Experimental  data  of  dosage  Fig.  4.  (1)  Rj,  R2?  (2)  Rq  diagram, 
distribution  of  gamma  radiation  during  Key:  (a)  m;  (b)  Kilotons;  (c)  1; 
detonation  of  Hiroshima-type  atomic  (d)  Ground  surface;  (e)  Air;  (f)  2. 
bomb:  1.  A*252  m;  2.  X«282  m;  3.  A= 

393  m. 


2.  Decayed  dosage  after  gamma  rays  pass  through  the  roof  of  an  air  defense 


basement 


According  to  the  practical  situation,  the  gamma  ray  incidence  acts  on  an 
oblique  direction.  First,  the  rays  pass  through  the  wall  of  a  story,  and  the 
story  itself;  then  the  rays  pass  through  the  cover  of  the  shelter-basement 
before  finally  entering  it. 

Calculation  formula:  IV*  *r Mh*!-**  *  (in 

From  Fig.  1,  we  obtain  a^»0.95;  from  Table  4,  we  obtain  (R'/H)*l  and  b*0.39. 

Converting  the  10-cm  covered  blast  furnace  slag  on  the  roof  of  shelter- 
basement  into  the  equivalent  thickness  of  reinforced  concrete 

« 1.1  cm. 

•  a  i.d 

The  wall  of  two-brick  thickness  is  converted  into  the  equivalent  thickness  of 
reinforced  concrete  IT»[(49x0.05)/0.07]*35  cm. 

Since  the  gamma  rays  are  of  oblique  incidence,  the  wall  thickness 

HnJ-(35/cos  45*)-49.5  cm. 

In  this  example,  a  six-story  building  is  above  the  shelter-basement;  it  is 
possible  that  gamma  rays  may  pass  through  five  floors  before  reaching  the 
basement  roof.  Therefore,  the  total  thickness  is 

Hw*l 0x5*3 . 9+30+49 .5*133.4  cm 

(Note:  the  thickness  of  a  story  floor  is  10  cm,  and  the  thickness  of  the  basement 
roof  is  30  cm.) 

Decayed  dosage:  d”*0.95x0.39x10700x2"13^‘*^10»0.38  roentgens. 

3.  Dosage  of  decayed  gamma  rays  after  passing  through  a  wall 

Calculation  for  mi  oTer,  arbDm2-/tcr*  <i*> 

Usually,  a  wall  oi  a.,  air  defense  shelter-basement  is  situated  below  the 
ground  surface;  only  a  wall  of  a  semi-basement  and  pile-up  type  air  defense  works 
is  exposed  to  the  ground  surface.  In  this  example,  this  portion  is  neglected. 

4.  Decay  of  gamma  rays  due  to  hole  opening 


a.  We  calculate  the  incident  radiation  dosage  according  to  the  theory 
of  angular  distribution 

As  revealed  by  experiments,  when  a  portion  of  the  radiation  dosage  is 
hindered  by  a  shielding  body,  the  incidence  of  rays  on  a  detector  is  only  the  gamma 
rays  exiting  from  a  scattered  medium  within  a  certain  three-dimensional  angle; 
this  leads  to  the  problem  of  angular  distribution  of  radiation.  It  follows  that 
the  incident  gamma  rays  onto  a  hole  opening  can  also  be  considered  as  incidence 
at  a  definite  three-dimensional  angle. 


Calculation  formula: 


(IS) 


As  recommended  by  related  manuscripts,  the  function  of  angular  distribution 
of  the  scattered  intensity  of  gamma  rays  in  the  air  due  to  fission  fragments 
can  be  written  as: 


src«>- r.'co.oosz  ♦  s.ia«'  *••«•»> 


(14) 


The  dosage  angular  distribution  of  the  captured  gamma  radiation  is  subject 


to 


F<i>  *  x,(o.oos  ♦  s.  it*1 


Cl.MN) 


(IS) 


F(0)  is  the  proportion  of  the  constituent  of  the  dosage  (in  the  total  dosage) 
in  a  unit  three-dimensional  angle  in  the  direction  of  the  scattered  angle  of 
the  constituent.  6  is  calculated  from  the  source  to  the  connecting  line  with 
the  detector. 


Kn  is  the  unitary  coefficient,  related  to  a  (the  intersection  angle  between 
the  ground  surface  and  the  direction  of  the  detonation  center)  and  earth's 
reflectivity.  *.■  ,  Jfcg- 

Kn  is  the  unitary  coefficient  of  the  assumed  black  medium  (for  earth)  without 
reflection  of  gamma  radiation. 

A  (<*)  is  the  integrated  reflectivity  of  the  earth  to  gamma  radiation  when 
Y° 

the  orientation  angle  of  the  detonation  point  is  a. 

A  major  portion  of  the  captured  gamma  radiation  has  been  reflected  within 

0.1  second;  its  function  actually  comes  to  an  end  within  0.3  second  following 
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detonation.  In  this  duration,  actually  the  radioactive  cloud  has  not  risen; 
therefore,  the  source  of  the  captured  gamma  radiation  can  be  considered  as 
a  stationary  source  like  the  source  of  neutrons. 

The  dose  rate  of  fragment  gamma  radiation  varies  with  time;  the  dose  rate 
is  determined  by  three  factors:  radioactive  decay  of  the  fragments,  rising  of 
radioactive  cloud,  and  a  vacuum  chamber  formed  by  the  shock  wave  in  the  atmo¬ 
sphere. 

The  dosage  of  gamma  radiation  within  a  unit  three-dimensional  angle  is: 


' ■«*,<!>*«<!** _ (lt) 

However,  6<u  -  arc cofon^rctt-^igj- 

-  a  ^aar>8co«*  ♦  cot  (arct|^pi^ 

-a)co«*j  (17; 

This  equation  is  too  cumbersome  in  application.  In  order  to  satisfy  the 
engineering  application,  the  author  proposes  to  use  the  following  formula. 

Calculation  formula ;  ,£  <u> 

In  the  equation,  is  the  .ddiation  dosage  (in  roentgens)  entering  the  hole 
opening;  when  ny-tel*  the  ratio  between  the  radiation  dosage  entering  the  hole 
opening  and  the  total  radiation  dosage  can  be  obtained  from  Fig.  5.  £2  is  the 
three-dimensional  angle  of  the  hole  opening  (£fca/r2);  a  is  the  cross  hatching 
of  the  three-dimensional  angle  of  the  hole  opening;  and  r  is  the  radius  (in 
centimeters)  of  the  three-dimensional  angle. 

For  the  example  given  before,  the  calculation  procedures  are  explained 
as  follows: 
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Due  to  rising  of  the  radioactive  cloud  following  detonation  of  an  atomic 
bomb,  at  any  instant  the  orientations  of  the  three-dimensional  angle  and  the 
detonation  center  vary;  i.e.,  the  intersection  angle  $,=$1+90®  is  the 
instantaneous  incident  angle) .  We  know  from  the  diagram  that  once  the  inter¬ 


section  angle  varies,  the  ratio  also  varies. 
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Fig.  5.  Ratios  between  the  radiation 
dosage  entering  the  hole  opening  and 
the  total  radiation  dosage  at  the  outer 
layer  of  the  works. 


According  to  conditions  for  this  example,  when 
radiation  is  0-20  percent,  $'»45°+90o*135°  and  ^*0.017; 
when  radiation  is  20-40  percent,  $'=135°  and  tu=0.017; 


when  radiation  is  40-60  percent,  $,=44o+90°=134°  and  nj=0.019; 


when  radiation  is  60-80  percent,  ^,=40°+90c=130®  and  n4=0.020; 


and  when  radiation  is  80-100  percent,  $,*32°+90°=122°  and  n5=0.022. 


0«  «/r  •  3  x  1*0  +  2  x  110  *  400*  •  0.058 
On  -  107*0  X  0.0M(  0.017  ♦  0.017  ♦  0.010 

♦  0.t*0*0.0ti)/8»  n.70  (roentgens) 


b.  The  incident  radiation  dosage  calculated  on  the  decay  coefficient 


The  gamma  rays  transmitted  to  the  underground  works  from  the  detonation 
center  move  along  straight  lines;  the  rays  should  be  turned  90°  through  the  hole 
opening,  thus  weakening  the  radiation  dosage.  There  are  differences  in  decay 
coefficients  in  different  papers;  as  recommended  by  the  INTRODUCTION  TO  NUCLEAR 
BOMB  DEFENSE  (published  in  the  United  States) ,  the  radiation  may  be  weakened  by 
0.07  fold;  in  some  papers,  it  is  suggested  that  the  radiation  may  be  weakened  by 


0.25  fold.  This  is  because  the  radiation  is  further  weakened  when  entering  into 
a  horizontal  passage  due  to  different  cross  sections  of  the  passage  and  the  inlet 
opening.  This  weakening  function  of  the  radiation  is  generally  called  as 
additional  function  or  area  function.  The  transmission  coefficient  AF  of  the 

ft 

area  function  can  be  expressed  by  the  following  equation. 

AFa*(area  of  horizontal  inlet  opening)/ (cross-sectional  area  of  horizontal 
passage)  (19) 

When  radiation  enters  the  underground  air  defense  structure  from  a  horizontal 
passage,  there  is  similarly  an  area  function,  as  expressed  by  the  following  formula 
AF^= (cross-sectional  area  of  horizontal  passage)/ (cross-sectional  area  of 
the  shelter  basement)  (20) 

It  is  apparent  that  radiation  intensity  is  weakened  from  passing  through  a 
protective  door.  As  sufficient  data  are  not  available,  the  author  proposes  for 
the  time  being  to  take  0.92  as  the  transmission  coefficient;  i.e.,  the  radiation 
intensity  Drj  of  gamma  radiation  after  passing  through  a  protective  door. 

Dy3=DYoxAFaxAFax^transiniss*on  coe^ici®nt  of  turn) x (transmission  coefficient 
of  the  door) 

According  to  the  conditions  of  this  example,  the  calculation  proceeds  as 
follows: 

Known  quantities:  cross-sectional  area  of  inlet,  90x160; 
cross-sectional  area  of  the  passage,  200x240;  and 

cross-sectional  area  of  the  shelter  basement,  240x600  (indicating  one 
room  of  the  shelter-basement) . 

i.  We  now  consider  that  there  are  three  90°  turns;  the  incident 

3 

dosage  is:  D^3*0.07  xl 0700*3.67  roentgens. 

ii.  AF  *(1.0x1 ,7)/ (2. 0x2. 4) *0,35. 

ft 

iii.  AF' *(2.0x2. 4)/ (2. 4x6.0) *0.33. 

a  17 


iv.  The  decayed  radiation  dosage  after  passing  through  10  cm  of 
reinforced  concrete 


Drt'.flTxDT,-2  *i»  =o.95*3.*r*2-,,'M  =1.75  roentgens. 


The  incident  dosage  of  the  hole  is: 

D  *D  xAF  xAF * x f  t ransmi s s i on  coefficient  of  turn)xD*'  =10700x0. 35x0. 33x 
y3  yo  _  a  a  Y3 

0.07+1.75*2.17  roentgens. 


v.  Dosage  of  gamma  radiation  entering  the  shelter -basement: 
D  =0.38+11.79+2.17=14.34  roentgens. 


IV.  Calculation  for  Defending  Against  Neutron  Stream  by  an  Air  Defense  Structure 

1.  Calculation  of  incident  dosage  DnQ  of  neutron  stream 

D  *(m/R2)(e“R/,25°)  (biological  roentgen-equivalent)  (/l) 
no 

In  the  equation,  m  -  is  the  decided  upon  coefficient  of 
A-bomb  equivalent  which  is  employed  according  to  Fig.  7. 

P„ » «***'“* -  12S0Q  (biological  roetgen  equivalent) 
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2.  Calculation  of  the  decayed  dosage  after  passing  through  the  roof  of 
shelter-basement  with  consideration  of  oblique  emission  of  neutron  stream 

Calculation  formula:  D,  =  aj)m2 co»,  (20 

In  the  equation,  is  the  dosage  after  a  neutron  stream  passes  through  the  roof 

of  shelter-basement  with  thickness  (of  roof)  of  H  X/cos  *  ;  d  is  the  half-value 

n  no 

thickness  of  the  neutron  stream,  obtained  from  a  table;  and  an  is  the  decay 
coefficient  for  neutrons  with  consideration  of  a  medium  surface,  obtained  from 
Table  8. 

a.  The  story  floor,  wall  and  blast-furnace  slag  are  converted  into 
equivalent  values  of  thickness  of  reinforced  concrete. 

For  10-cm  blast  furnace  slag:  H^#=(10x9.7)/(11,6)=8.4  cm; 

for  wall:  ^=(49x12)/  (11.6)*50.7  cm;  and 

for  story  floor;  H'" *5x10=50  cm. 
na 

The  roof  of  shelter-basement  is  30  cm  thick.  Total  thickness  (of  cover): 

H  =8 . 4+50. 7+50+30=139  cm. 
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Fig.  7.  q<Mn  diagram. 
Key:  (*)  Kilotons. 


Table  8. 


(c)  ±*. 


(e)«  « 


Key:  (a)  Material;  (b)  Thickness  of  half¬ 
value  decay;  (c)  Soil,  dressed  masonry;  (d) 
Concrete;  (e)  Iron  or  steel. 
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c.  Decayed  dosage:  I>M  =  0.8x12500  xj-“»'*’‘*;I-».o.062 

(biological  rpentgen-equivalent) 

2.  Calculation  of  decayed  dosage  after  passing  through  a  wall  with  considera 
tion  of  oblique  emission  of  neutron  stream: 


CTO 


Calculation  formula: 


£>„  =  a,Dm  2‘*‘*  ilt»‘ 


<ii> 


Since  the  shelter-basement  is  situated  below  the  ground  surface,  this  portion 
of  incident  dosage  is  neglected. 

3.  Decayed  dosage  of  neutron  stream  entering  a  hole  opening 

a.  Dosage  of  neutron  stream  entering  through  a  door  opening  calculated 
with  consideration  of  angular  distribution: 

Calculation  formula:  (22) 

(To  be  continued  on  back  cover  of  text) 

[Figure  6  and  Table  7  are  missing;  they  may  be  in  the  last  portion  of  text  that 
was  missing.  Translator's  note] 


(Continuation  from  text  page  32) 

It  is  now  assumed  that  the  vertex  of  the  three-dimensional  angle  is  located 
at  the  midpoint  of  two  passages,  assuming  that  r=120  cm  and  the  included  angle 
*’=135°.  From  Fig.  5,  nn=0.012,  s=2xl00+2xl80=580  cm,  and  n=(580/1202)=0.04. 

D^3=12500x0. 012x0. 04=6  (biological  roentgen-equivalent) 

(b)  Dosage  of  decayed  neutron  stream  after  passing  through  the  pro¬ 
tective  door 

/ 

Calculation  formula:  D  =a  D',.2-h/dn  =  0.6x6.0*2"10/11 ,6  =1.98 

n2  n  n3 

(biological  roentgen-equivalent) 

4.  Total  dosage  of  neutron  stream  entering  the  basement  shelter 

Dn=0. 062+0+1. 98=2. 042  (biological  roentgen-equivalent) 

5.  As  derived  from  the  aforementioned  steps,  it  is  required  that  the 
summation  of  the  remaining  total  dosage  of  gamma  rays  and  neutron  stream  should 
not  exceed  the  allowable  dosage,  50  roentgens. 

That  is,  D=DY+Dn=14. 34+2. 042=16. 382  [50] 

V.  Characteristics  in  Calculating  Penetrating  Radiation  After  Detonation  of 
a  Thermonuclear  Weapon 

As  mentioned  above,  the  existence  of  ultrafast  neutrons  is  a  characteristic 
>  of  penetrating  radiation  upon  detonation  of  a  thermonuclear  weapon.  For  various 
types  of  hydrogen  bombs,  gamma  rays  still  occupy  the  major  portion  of  the 
radiation;  generally,  fast  neutrons  and  ultrafast  neutrons  do  not  occupy  a  major 
portion.  The  following  formulas  are  used  in  calculation: 

13  2  -R/250 

Dosage  of  gamma  rays:  Dy*(4xl0  /R  )*e  (roentgens)  (23) 

12  2  -R/250 

Dosage  of  fast  neutrons:  Dn=(3xl0  /R  )*e  (biological  roentgen- 

equivalent)  (24) 

Dosage  of  ultrafast  neutrons:  D  jtra>  (4xl011/R2)  *e  R<^°  (biological  roentgen- 

equivalent)  (25) 
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VI.  Conclusion 


Through  the  above  discussion,  the  authors  have  formed  an  understanding 
of  calculation  of  penetrating  radiation  defense  in  fortifications.  The  following 
points  are  listed  for  the  readers'  attention  to  avoid  mistakes. 


1.  In  this  paper,  it  is  assumed  that  the  ascending  speed  (of  the  radio¬ 
active  cloud  following  a  nuclear  detonation)  is  constant;  that  is,  Ht«H+v  t 

(v  =100  m/sec).  As  stated  by  associate  professor  Qian  Shihu,  who  quoted  a  book, 
cp 

Protection  Against  Neutron  Stream  and  Gamma  Radiation  in  Nuclear  Radiation, 

1969,  written  by  Kuhetaiweiji  [transliteration  of  Russian  name]  et  al  of  the 

Soviet  Union  that  the  ascending  radioactive  cloud  has  constant  acceleration; 

2 

i.e.,  Ht=H+6t  .  Since  there  is  no  such  statement  in  the  recently  published 
literature  in  the  United  States  and  Britain,  the  constant  acceleration  of  the 
ascending  radioactive  cloud  is  not  included  (authors'  note). 


2.  In  this  paper,  the  ascending  radioactive  cloud  has  constant  velocity; 
therefore,  the  incident  dosage  is  calculated  by  using  Eq.  (3).  If  the  ascending 
cloud  has  constant  acceleration,  then  the  incident  dosage  is  a  variable.  The 
value  of  Dyo  ^  can  be  derived  by  integration.  Any  measured  data 

can  be  used  directly. 


3.  When  the  consideration  is  given  that  the  ascending  radioactive  cloud  has 
constant  acceleration,  data  in  Tables  4  and  5  can  still  be  used  as  an  approxima¬ 
tion,  the  authors  will  compile  blank  forms  of  the  oblique-direction  decay 
coefficient  of  the  protective  layer. 

As  the  authors  have  small  scientific  learning  with  limited  level,  mistakes 
naturally  exist  in  the  paper.  The  authors  welcome  corrections  from  experts, 
scholars  and  readers. 
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